Spin torque from tunneling through impurities in a magnetic tunnel junction 
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We calculate the contribution to the spin transfer torque from sequential tunneling through im- 
purities in a magnetic tunnel junction. For a junction with weakly polarized ferromagnetic contacts, 
the torque is found to be in the plane spanned by the magnetizations of the ferromagnetic contacts 
and proportional to sin#, where 9 is the angle between the magnetic moments. If the polarization 
is larger, the torque acquires a significant out-of-plane component and a different dependence on 9. 



PACS numbers: 85.75.-d, 72.25.-b, 73.23.Hk 



I. INTRODUCTION 



A magnetic tunnel junction consists of two ferromag- 
netic layers, separated by an insulator.— i 2 - Magnetic tun- 
nel junctions play a prominent role in proposals for 
magnetic memory applications, where the information 
is stored in the relative orientation of the magnetiza- 
tions of the two ferromagnetic layers. The two essen- 
tial operations for the realization of a memory element, 
reading and writing, rely on the tunneling magnetoresis- 
tance effect (the observation that the resistance of the 
junction depends on the relative orientation of the two 
ferromagnets 3 -) and the spin transfer torque (the fact that 
passing a large current through a magnetic tunnel junc- 
tion exerts a torque on the magnetizations of the two 
layers4^) Although both effects also exist in metallic 
magnetic junctions, tunneling junctions have a signifi- 
cantly higher magnetoresistance and an impedance that 
is better matched to the requirements of semiconductor 
technology.— 

Although the tunneling magnetoresistance effect has 
been studied experimentally and theoretically for more 
than 30 years, the spin transfer effect in magnetic tunnel 
junctions had not received attention until recently. The 
spin transfer torque consists of two contributions, a com- 
ponent in the plane spanned by the magnetization direc- 
tions of the two ferromagnetic layers in the junction, and 
a component transverse to that plane. The out-of-plane 
component is also referred to as the "field-like" torque 
or the "nonequilibrium exchange interaction" . The in- 
plane torque component can be understood in terms of 
the non-conservation of spin for electrons tunneling be- 
tween the ferromagnetic layers^ The out-of-plane torque 
component has been attributed to a more subtle in- 
terference between electrons reflected off the front and 
back ends of the insulating barrier*^ Theoretical pre- 
dictions exist for both torque contributions, with similar 
conclusions based on tight-binding calculations that in- 
clude the full band structure of the junctio n 10 ' 11 : 12 and 
free electron models that neglect most band structure 
cffcctsj£&i£ On the experimental side, quantitative mea- 
surements of both components of the spin transfer torque 
in tunnel junctions are possible using spin-transfer-driven 
ferromagnetic resonance i 14 ' 15 



Whereas most of the theoretical approaches to the 
spin-transfer torque in magnetic tunnel junctions ad- 
dress ideal junctions, impurities and defects are known 
to be abundant in the MgO magnetic tunnel junc- 
tions used in experiments ! 16 ' 17 : 18 While impurities in- 
side the barrier are known, both theoreticall y 19 ' 20 and 
experimentally, 18,21 to strongly affect tunneling magne- 
toresistance, there has been less attention to their effect 
on the spin transfer torque . 13 ' 22 

Unlike for a theory of spin torque and magnetocon- 
ductance in ideal junctions, where both phenomena are 
attributed to electron states that are extended on both 
sides of the junction, impurity-mediated transport in- 
volves localized states inside the barrier. Inside the bar- 
rier, electron-electron interactions are poorly screened, 
and interactions cannot simply be accounted for on the 
mean-field level. In this letter, we describe the impurity- 
mediated transport using a model of sequential tunnel- 
ing that is able to incorporate interaction effects exactly. 
The sequential tunneling picture is valid for high tem- 
peratures or voltages (temperature and/or voltage larger 
than the width of the impurity state), which makes our 
theory applicable to all but the thinnest junctions. This 
puts our theory in an entirely different parameter range 
than Refs. and[H, which addressed zero temperature 
torques in the absence of interactions. 

The concept of sequential tunneling has been used 
previously to describe spin-dependent transport through 
metal nanoparticles (or "quantum dots" ) tunnel coupled 
to ferromagnetic electrodes! 23 ' 24 ' 25 ' 26 For these systems, 
one considers how a current between the ferromagnetic 
contacts affects the spin accumulated on the nanopar- 
ticle. In the present context, we are interested in the 
opposite effect, the back-action of the current on the 
polarization of the ferromagnetic contacts. Below, we 
first analyze the current I and torques T^, and on 
the two ferromagnets in the tunnel junction for sequen- 
tial tunneling through a single impurity with a single 
(spin degenerate) energy level, following ideas laid out 
in the theory of spin-dependent transport through quan- 
tum dots. We then turn to the case of magnetic tunnel 
junctions, where impurity-mediated contributions to the 
current and torque have to be summed over many impu- 
rities in the junction. 
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II. SEQUENTIAL TUNNELING THROUGH A 
SINGLE IMPURITY 

We denote the energy of the singly-occupied and 
doubly-occupied impurity level by e and 2e + U , re- 
spectively, where the interaction energy U accounts for 
the Coulomb repulsion of the two electrons in a doubly- 
occupied impurity site. The left (L) and right (R) fer- 
romagnets are held at chemical potential /i£ and [i Rl 
respectively. For definitencss, we assume that the bias 
eV = Hl — Mi? > 0. (See Fig. Q] for a schematic pic- 
ture.) The direction of the magnetization of each fer- 
romagnetic contact is denoted by the unit vector m Q , 
a = L, R. Tunneling between the impurity level and the 
reservoirs is described by tunneling rates Tlct and T RlT , 
where a =|, J, for the majority and minority spin direc- 
tions, respectively. 

The theory of sequential tunneling is applicable if the 
temperature T is much larger than the tunneling rates 
T aa . In this regime, interference effects (such as those 
responsible for an impurity-mediated zero-bias torque 2 ^) 
are smeared out and impurity-mediated transport can 
be described using the probabilities pj that the impu- 
rity level is occupied by j electrons, j — 0, 1, 2, and the 
expectation value s of the impurity spin. 26 Their time 
derivatives pj and s are expressed in terms of particle 
and spin currents onto the impurity site, 

Po = -I'l-I'r, 

Pi = I'l+I'r-I'L-Ir, s = I L +I R , 

p-2 = il + i'k, (i) 



where I' a and 1'^ are the currents from contact a for pro- 
cesses in which the occupation changes from to 1 and 
from 1 to 2, respectively, and I Q is the spin current onto 
the impurity site, a — L, R. The component of I a per- 
pendicular to m Q causes a torque Tq4, 



T a = h(I a x m Q ) x m„ . (2) 



The torques are decomposed into an in-plane component 
T a n, a = L, R, and an out-of-plane component Tj_, 

T L = T L |,m L x n + T ± n, (3) 
T^ = T m m R x n -Tj.il, (4) 

where n is the unit vector perpendicular to and m/j, 



n = (m L x m R )/\m L x m R \. (5) 



Following Ref. [H, the currents I' a , 1'^, and l a are ex- 
pressed in terms of the tunneling rates T aa and the dis- 
tribution functions f a (0 = [1 + exp((£ — /x Q )/T)] _1 of 
the two reservoirs, 



I' a = 2T a+Po f a (e)- [1 -/„(£)] (r a+ pi + 2r Q _m tt -s), 

I'L = fa(e + U) (T a+ - 2r Q _m Q • s) - 2T a+ p 2 [l - f a (e + (7)], 

I Q = h a m a x s + r a -f a (e)p m a - T a -[1 - f a (e + U)}p 2 m a 

- [1 - /«(£)] (r Q+ s + T a ^ Pl m a /2) - f a (e + U) {T a+ s - r Q _ Pl m Q /2) . 



Here r a± = (r aT ± T a {)/2, and 



2tt 



in 



j7(i-2/ a (Q)r a _ 

[e-Z)(e + U-Z) 

£ - Ma 



e + U - [i a 



(6) 



is an effective exchange field arising from virtual tun- 
neling from the impurity level to the reservoirs^ The 
exchange field is an interaction effect: h a = in the 
absence of Coulomb repulsion on the impurity site. 

Sequential tunneling takes place if at least one of the 
energies e and e + U lies between \ih and /i R , see Fig. 
[TJ Since typical interaction energies U are large (up to 



several eV), we assume that eV <C U. In this regime, the 
occupation of the impurity site can change by at most one 
electron. Depending on whether e or e + U lies between 
the chemical potentials of the source and drain contacts, 
the occupation of the impurity site fluctuations between 
and 1 ("case I") or 1 and 2 ("case II"). In case I, the 
only nonzero charge currents are I' L and I' R , whereas in 
case II, the nonzero charge currents are I'l and 1'^. Be- 
low, we give the relevant expressions for case I only. Case 
II follows from the expressions below by making the sub- 
stitutions L <-> R and h «-» —h. 

Solving Eqs. ©-© for T < eV, we then calculate 
the current I = II through the impurity as well as the 
torques T^ and T R . We then find 
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/ = D~ l Y L+ [(r| + +hl + h% + 2h L h R cos9)(T 2 R _ - T% + ) - r 2 R _h 2 L sm 2 6] , (7) 
T m = 2hD~ l sin0 {T L .{Y 2 R _-T 2 R+ )[Y 2 R+ + {h 2 R + h\ + h L h R co*9)\- h 2 L Y R _(Y R _Y L _ -Y L +Y R +cos9)} , (8) 

Tin - 2n J D- 1 /i L sin0[r L _r2j_(^ + ^LCos0)-r fi+ (/ lL r L+ r fl _ + /i fl r fl+ r L _)] , (9) 

T x = 2hD- 1 h L smO {T L -T R+ h L h R -Y L+ Y R ^(Y 2 R+ + h 2 R ) - [Y L+ Y R _h L h R -T L _T R+ (T 2 R _ + h 2 R )\cosQ} , (10) 
with 

D = 2{Y 2 R+ + h 2 L + 2h L h R cos9 + h 2 R )[Y R _(Y R _+ 2T L _cose)-r R+ (r R+ + 2T L+ )} 

+ 2T R -h L {2r L -h R -r R -h L )sm 2 8. (11) 



It is important to point out that the in-plane torque 
components Tr|| and Tn\ on the magnetizations of the 
source and drain contacts are, in general, not equal. (The 
out-of-plane components are equal because of conserva- 
tion of angular momentum.) The origin of this effect 
is the Coulomb interaction on the impurity site, which 
lifts the symmetry between source and drain reservoirs. 
To see this, recall that in case I (for which the above 
equations are derived), the Coulomb interacation forbids 
double occupancy of the impurity site. Since the electron 
spin is not changed when an electron tunnels onto an 
empty impurity level, the electron spin may be changed 
only when the electron tunnels off the impurity site, or 
because of the action of the exchange fields Hl and h R . 
This explains why T R \\ is nonzero to lowest order in the 
exchange fields, whereas T L \\ and Tj_ are of higher order 
in }il and h R , respectively. 

III. IMPURITY-MEDIATED TUNNELING IN A 
MAGNETIC TUNNEL JUNCTION 

In a magnetic tunnel junction there will be many im- 
purities, distributed spatially throughout the insulating 
spacer layer, and with a distribution of energy levels s. 
(The energy e depends on the electrostatic environment 
of the impurity, hence its fluctuations.) In order to add 
the contributions from all impurities, we characterize the 
impurity configuration by the density p{e 1 x) of impuri- 
ties per volume and per energy, where < x < d is the 
distance from the source reservoir, d being the width of 
the insulating spacer layer (see Fig.Q]). We neglect fluctu- 
ations of the interaction energy U, which is less suscep- 
tible to the electrostatic environment of the impurities 
than e. 

We assume that the density of impurities is small 
enough, 

hpd 3 T < 1, (12) 

where hY is the typical width of an impurity energy 
level, so that the transport of an electron will take place 
through a single impurity only, i.e., multiple-impurity 
processes are ignored. The total current I and torque T 



are then written as integrals over the energy e of impu- 
rities in the junction, 

/= f R de[t'(e)+t"(e)], T= f * ds[t'(s) +t"(e)}. 
J Jul 

where i', t' and i" , t" represent contributions from impu- 
rities at energy e (case I) and e — U (case II), respectively. 

For barriers with rough interfaces the transverse mo- 
mentum is not be conserved at the ferromagnet-insulator 
interface. In that case, the tunneling rates Y aa can be 
estimated as£ 

Fqct Pair & ^ j (1^) 

where P aa is a polarization factor and A the wavefunc- 
tion decay length in the spacer layer 3& We will take the 
junction to be symmetric, Pj, CT = P Ra — P a , and define 
P T = (Pj =F -Pj.)/2- For barrier width d 3> A the current 
and torque are dominated by impurities near the cen- 
ter of the barrier, x « d/2J£ The contribution of these 
impurities to the current and torque scales ~ e~ rf / A if 
d 3> A, whereas there is a faster exponential suppression 
for direct tunneling or for impurities located off-center. 

We have calculated the two contributions to the cur- 
rent and torque for arbitrary angle and polarization 
p = P-/P+. The results are rather lengthy, and we re- 
fer to Fig. [3 for a numerical evaluation for a few repre- 
sentative values of p. Closed form expressions could be 
obtained in the limit of weakly polarized ferromagnets, 
p < 1 only. For impurities with p R < e < (Xl (case I), 
we find 

i' = i' (l + ^p 2 cos 6 + ^p 5/2 vl sin 2 9 + . . . J , 
t' m = t / sm9^l-^p 1 / 2 r lL cos9 + , 
t' L || = sin 9 (2P 1 ' 2 - p 3 / 2 7f R cos 9 + . . .) , 

t'x = |Vsin#(2p 1 / 2 +p 3 / 2 f^cos0 + ...) , (15) 
where 

i' Q (e) = ir2- 7 / 2 P + \p(s, d/2)e- d /\ t' (e) = 2hpi' (e), 

(16) 
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and rj a = y/\h a \/T a -. In Eq. ([T5|l we kept only those 
sub-leading terms in the small-p expansion that have a 
different 9 dependence than the leading terms. The case 
II contributions are obtained by interchanging L and R 
and replacing p(e, d/2) by p(e — U, d/2). The total torque 
is found by adding contributions for cases I and II and 
integrating over e. 

There is a striking asymmetry between the impurity- 
mediated in-plane torques on the source and drain reser- 
voirs. For small polarization (jp <C 1) and case I, the 
in-plane torque tm dominates over the other two torque 
terms. This is markedly different from the torque from 
direct tunneling, which has equal magnitudes for to and 
tR\\Z The situation is reversed for case II. However, as 
there is no a priori reason why the spectral impurity den- 
sities at energies e and e — U, which set the sizes of the 
torques for cases I and II, are equal, one still expects 
the magnitudes to and tm to be rather different after 
torques from the two cases are added. 

It is important to point out that both the source- 
drain asymmetry and the existence of an out-of-plane 
component of the torque are interaction effects. The 
source-drain asymmetry was already discussed at the end 
of Sec. [ill That the out-of-plane torque is an interac- 
tion effect can be seen explicitly from Eq. ([TO)) , which 
is proportional to {Hr for case II). Without inter- 
actions, h L — h R = 0. (Note that the out-of-plane 
torque in theories without electron-electron interactions 
has an altogether different origin: It is caused by inter- 
ference effects For the temperature range we con- 
sider, T 3> KT, interference effects are smeared out for 
impurity-mediated transport.) 

For strongly polarized ferromagnetic contacts (p ~ 1) 
all three torque terms are of comparable magnitude, al- 
though to and tfjy remain different. The order of magni- 
tude of the torque, normalized to the impurity-mediated 
current i' , is the same as for the case of direct tunneling, 
normalized to the direct tunneling current 7 -. Whereas 
the angular dependence of all three torque terms is the 
standard geometric sin 9 dependence characteristic of the 
direct tunneling torque 7 -^! 9 - ' 11 ! 12 if p <C 1, the 9 depen- 
dence is more complicated for strongly polarized contacts 
because of the presence of the exchange fields and Iir, 
see Fig. 



IV. DISCUSSION 

The impurity-mediated torque considered here coex- 
ists with the torque from direct tunneling. How can 
the two contributions to the torque be separated ex- 
perimentally? The calculation of the previous sections 
shows that there are four significant differences between 
the impurity-mediated torque and the direct-tunneling 
torque: (i) The impurity-mediated torque is not sym- 
metric under reversal of the bias, whereas the torque 
from direct tunneling is symmetric under bias reversal. 
The same holds for the impurity-mediated and direct- 



tunneling currents. Hence, the presence of a bias asym- 
metry for the current or torque is an indicator for the 
order-of-magnitude of the impurity-mediated contribu- 
tion, (ii) The impurity-mediated torque has a non- 
trivial angular dependence for strongly polarized con- 
tacts, whereas the direct-tunneling torque has a sinu- 
soidal angular dependence. The angular dependence of 
the torque is observable in current-induced ferromagnetic 
resonance experiments ! 14 > 15 (iii) Impurity-mediated and 
direct-tunneling contributions to the current and torque 
have different dependences on the thickness d of the insu- 
lating layer: The impurity-mediated contribution scales 
oc e~ rf / A , whereas the direct-tunneling contribution scales 
cx e~ 2d / x . This means that the impurity-mediated torque 
will dominate for sufficiently thick junctions, irrespective 
of the impurity concentration, (iv) The bias-dependence 
of the impurity-mediated torque tracks the spectral im- 
purity density p(e, x), whereas the bias-dependence of the 
direct-tunneling torque should be less pronounced. (The 
spectral impurity density is not the only source of a bias 
dependence: One also expects a bias dependence from 
the bias-dependence of the tunneling rates T. However, 
that bias dependence is likely to be the shared between 
the two contributions to the torque.) For the impurity- 
mediated torque, resonant features in the spectral den- 
sity, which can be brought out using more detailed the- 
oretical modelin g 19 i 20 ' 22 or additional experiments ; 18 ^ 1 
can then be correlated with the bias dependence of the 
torque. 

It is instructive to give a rough comparison of the mag- 
nitudes of the impurity-mediated and direct tunneling 
currents or torques. Since they have different exponen- 
tial dependences on the barrier thickness d, we estimate 
the pre-exponential factors by order of magnitude only 
and neglect differences between the three relevant micro- 
scopic length scales (Fermi wavelengths for majority and 
minority electrons in the ferromagnetic contacts, wave- 
function decay length A). Dimensional analysis then esti- 
mates the ratio of impurity-mediated and direct currents 

^impurity: -^direct and torques Tj m p Ur ity, ^direct ^ 

-^impurity /-^direct ~ ^impurity /^direct ^ 6 ^ TlX , (17) 

where n = J depie) is the impurity concentration. Ref- 
erence uses n ~ 10 27 m~ 3 , whereas an estimate using 
the barrier height in MgO and the effective mass of the 
electron gives^ 7 - A ~ 10 -9 m. From this, we conclude 
that impurity-mediated transport dominates already for 
barrier thickness d > 10 -9 m, well within the experi- 
mentally relevant range. A more accurate comparison 
requires knowledge about the spectral density specific to 
the impurity type and is beyond the scope of this article. 

Acknowledgments 

We acknowledge helpful discussions with Bob 
Buhrman, Dan Ralph, John Read, Joern Kupferschmidt 
and Lin Xue. This work was supported by the Cornell 



5 



Center for Nanoscale Systems under NSF Grant No. 0705476. 
EEC-01 17770 and by the NSF under Grant No. DMR 



X.-G. Zhang and W. H. Butler, J. Phys.: Condens. Matter 
15, R1603 (2003). 

E. Y. Tsymbal, O. N. Mryasov, and P. R. LeClair, J. Phys.: 

Condens. Matter 15, R109 (2003). 

M. Julliere, Phys. Lett. A 54, 225 (1975). 

J. C. Slonczewski, J. Magn. Magn. Mater. 159, 1 (1996). 

L. Berger, Phys. Rev. B 54, 9353 (1996). 

S. S. P. Parkin, K. P. Roche, M. G. Samant et al., P. M. 

Rice, R. B. Beyers, R. E. Scheuerlein, E. J. O'Sullivan, 

S. L. Brown, J. Bucchigano, D. W. Abraham, et al., J. 

Appl. Phys. 85, 5828 (1999). 

J. C. Slonczewski, Phys. Rev. B 71, 024411 (2005). 

J. Xiao, G. E. W. Bauer, and A. Brataas, Phys. Rev. B 

77, 224419 (2008). 

A. Manchon, N. Ryzhanova, A. Vedyayev, M. Chschiev, 
and B. Dieny, J. Phys.: Condens. Matter 20, 145208 
(2008). 

C. Heiliger and M. D. Stiles, Phys. Rev. Lett. 100, 186805 
(2008). 

I. Theodonis, N. Kioussis, A. Kalitsov, M. Chshiev, and 

W. H. Butler, Phys. Rev. Lett. 97, 237205 (2006). 

A. Kalitsov, I. Theodonis, N. Kioussis, M. Chshiev, W. H. 

Butler, and A. Vedyayev, J. Appl. Phys. 99, 08G501 

(2006). 

A. Manchon, N. Ryzhanova, N. Strelkov, A. Vedyayev, and 

B. Dieny, J. of Phys.: Condens. Matter 19, 165212 (2007). 

A. A. Tulapurkar, Y. Suzuki, A. Fukushima, H. Kub- 
ota, H. Maehara, K. Tsunekawa, D. D. Djayaprawira, 
N. Watanabe, and S. Yuasa, Nature 438, 339 (2005); 
H. Kubota, A. Fukushima, K. Yakushiji, T. Naga- 
hama, S. Yuasa, K. Ando, H. Maehara, Y. Nagamine, 
K. Tsunekawa, D. D. Djayaprawira, et al., Nature Physics 
4, 37 (2008). 

J. C. Sankey, Y.-T. Cui, J. Z. Sun, J. C. Slonczewski, R. A. 
Buhrman, and D. C. Ralph, Nature Physics 4, 67 (2008). 
P. G. Mather, J. C. Read, and R. A. Buhrman, Phys. Rev. 
B 73, 205412 (2006). 

J. J. Cha, J. C. Read, R. A. Buhrman, and D. A. Muller, 
Appl. Phys. Lett. 91, 062516 (2007). 
G. X. Miao, Y. J. Park, J. S. Moodera, M. Seibt, G. Eilers, 
and M. Miinzenberg, G. X. Miao et al., Phys. Rev. Lett. 
100, 246803 (2008). 

E. Y. Tsymbal, A. Sokolov, I. F. Sabirianov, and 

B. Doudin, Phys. Rev. Lett. 90, 186602 (2003). 

J. P. Velev, K. D. Belashchenko, S. S. Jaswal, and E. Y. 

Tsymbal, Appl. Phys. Lett. 90, 072502 (2007). 

R. Jansen and J. S. Moodera, Appl. Phys. Lett. 75, 400 

(1999); Phys. Rev. B 61, 9047 (2000). 

M. Y. Zhuravlev, E. Y. Tsymbal, and A. V. Vedyayev, 

Phys. Rev. Lett. 94, 026806 (2005). 

X. Waintal and P. W. Brouwer, Phys. Rev. Lett. 91, 
247201 (2003). 

A. Cottet, W. Belzig, and C. Bruder, Phys. Rev. Lett. 92, 
206801 (2004). 

J. Konig and J. Martinek, Phys. Rev. Lett. 90, 166602 
(2003). 



Konig, and J. Martinek, Phys. Rev. B 70, 




FIG. 1: Schematic picture showing the impurity energy level 
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The left panel is for the case (j,r < e < [il, in which the 
occupation of the impurity site fluctuates between and 1. 
The right panel is for the case /j,r < e + U < ^l, in which the 
impurity site occupation varies between 1 and 2. These two 
cases are labeled I and II in the main text. 
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FIG. 2: Spectral current i'(s) and torque components t' R ^(s), 
<4||(e), and t' ± (e) versus the angle 6 between the magnetiza- 
tions of the left and right ferromagnetic contacts. The curves 
shown are for polarization p — 0.2 (solid), 0.5 (short dash), 
0.9 (dash-dot), and 0.99 (long dash). For definiteness, the 
exchange field parameters have been set to t]l = tjr = 1. 
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